Construction of a dispersion relation from an arbitrary density of states 
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The dispersion relations of energy bands in solids are characterized by their density of 
states, but a given density of states may originate from various band structures. We show 
how a spherically symmetric dispersion can be constructed for any one-band density of 
states. This method is applied to one-, two- and three-dimensional systems. It also serves 
to establish that any one-band spectrum with finite bandwidth can be obtained from a 
properly scaled dispersion relation in the limit of infinite dimensions. 
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1. Introduction 

The band structure of electrons describes their quantum-mechanical motion in the 
periodic potential of a crystal lattice in the absence of interactions. This picture 
is modified by the Coulomb interaction between them, which must be taken into 
account accurately to understand the phenomena observed in strongly correlated 
systems such as transition-metal oxides, ferromagnetic metals, or the cuprate super- 
conductors. Nonetheless the properties of the non-interacting lattice system remain 
an important input quantity. 

The electronic spectrum of a solid is characterized by the density of states, 



where the integral is over the first Brillouin zone, which has volume Vb', here and 
in the following we consider only a single separated band eg. The mapping from 
dispersion relation to density of states is not invertible, i.e., two different dispersions 
may lead to the same function N(u>). Nevertheless in certain situations it is desirable 
to construct a dispersion that corresponds to a particular N(lo). In this paper we 
present a method that performs this inversion in any dimension D, with additional 
assumptions on the functional form of eg. 

Construction of a dispersion with given spectral properties is especially of inter- 
est within the framework of dynamical mean-field theory (DMFT) for Hubbard-type 




(1) 



models, which becomes exact in the limit of infinite dimensions, D 



oo.l 




In this 
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approach typical tight-binding dispersions, arising from hopping of electrons be- 
tween orbitals on neighboring lattice sites, lead to a density of states which extends 
to infinite values o£ u>. For example, on the hypercubic lattice N(u>) becomes a 
Gaussian bell curves 



iV H = 7^ ex p(-— )■ ( 2 ) 



Exponential tails are also obtained for the generalized face-centered cubic lattices 
and the generalized diamond lattice.u 

The DMFT is particularly well suited to describe the Mott-Hubbard metal- 
insulator transition.Q In order to work with a finite bandwidth, which better models 
finite-dimensional systems, DMFT studies of this transition (for recent work see 
Ref. ^| and references therein) have concentrated on next-neighbor hopping on the 
Bethe lattice in the limit of infinite connectivity. It leads to a semi-circular density 
of statesE 



N(u) = -^JI-uj 2 , |w|<l, (3) 
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which in this case is defined via the one-particle Green's function, since there is 
no underlying periodic crystal lattice and hence no dispersion relation. The Bethe 
lattice has hitherto been the only known infinite-dimensional system with a fi- 
nite bandwidth. Below we will show how a dispersion can be constructed on a 
Z?-dimensional crystal lattice that in the limit D — > oo yields any given density of 
states with finite bandwidth, providing an answer to one of the open conceptual 
questions about this limit. 

In DMFT all one-particle quantities in the homogeneous phase involve the dis- 
persion only implicitly through N(w). Therefore another standard approach has 
been to use a function N(u>) that derives from a dispersion on a finite-dimensional 
lattice or is tuned by hand. For instance, Wahle et alB introduced the model density 
of states, 



. , 1 + \/l - a 2 \f\ - up- . , 



to study the effect of asymmetry in N(u>) on the stability of ferromagnetism in 
the Hubbard model. Here a 6 [— 1; 1] is a tunable parameter; (^) reduces to the 
symmetric function (||) for a = 0. For a > spectral weight is shifted towards the 
lower band edge, while for a = 1 N(lu) has a square-root singularity at w = — 1 
similar to the tight-binding density of states for the fee lattice. Below we will derive 
dispersions in D = 1, 2, 3, oo that represent this density of states. 
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2. Isotropic Fermi surfaces on Bravais lattices 

We construct a dispersion with a given density of states as follows. In terms of the 
latter the band filling n{ji) for Fermi energy \x is given by 

/< 

n(/i) = / du N(v) , (5) 



which lies between and 1. Suppose now that the dispersion relation is spherically 
symmetric, = e(|fc|), with \k\ 2 = ki + • — hkp. Furthermore we assume that e(fc) 
increases with k; then the Fermi sea is simply connected and includes the point k 
= 0, and there is a one-to-one relation between Fermi energy fi and Fermi vector k. 
The band filling can thus be obtained in another, geometric way as the normalized 
volume of the Fermi sea, Vo(k), defined by 

k 

v D {k)= fdk's D (k'), (6) 



o 



where sc(fc) is its normalized surface, 



^JdkvJdk D SUe-^ 



s D {k) = — dkv dk D 5[k-Jk 2 + --- + k 2 D ). (7) 



The dispersion e(k) is now defined by setting vo{k) equal to the band filling calcu- 
lated from the density of states, i.e., 

v D (k) =n(e(k)). (8) 

As is appropriate for a one-band density of states we assume N(lu) > for oj m i n < 
u> < cij max , so that n(p) is invertible on this interval. In terms of its inverse function 
/i(n) we can write 

e(k) = fi(v D (k)) . (9) 

Note that the functional form of n(n) depends only on the given density of states 
N(ui), whereas that of vo(k) depends only on the crystal lattice. Once the latter 
has been calculated for a lattice of interest, e(k) is immediately available for any 
N(to). 

Similarly, the density of states can be at once calculated from a given spherically 
symmetric dispersion e(fc), provided t'{k) > 0. In terms of fc(e), the inverse function 
of e(fc), we have n(e) = vjj(k(e)) from (^), and hence 

N(e) = n'(e) = s D (k(e))k'(e). (10) 



To apply the method in either direction it remains to calculate the functions su(fc) 
and vn{k). 
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3. Dimensions D = 1, 2, 3 

We consider now the special case of a hypercubic lattice in D dimensions, treating 
explicitly the cases D = 1,2,3, as well as the limit D — > oo. For convenience we 
measure A; in units of 2n/a, where a is the lattice spacing, so that the first Brillouin 
zone is located in the region |jfej| < 4, and the maximum value for \k\ is y/D/2. 
For k < | the integrations in (]?]) can be performed in spherical coordinates, since 
the vector k always remains within the first Brillouin zone. In this case we have 
s D (k) = Sok - 1 and v D {k) = V D k D , where S D = 2tt d / 2 /T{D /2) and V D = 
■k d I 2 /r(l + D/2) are the surface and volume of the D-dimensional unit sphere, 
respectively. This case applies for all k in D = 1, so that 

si(fc)=2, Vl (k)=2k, for Jfe e [0; §] , (11) 

representing the contribution of the two Fermi points. The spherical symmetry 
reduces to inversion symmetry in this simple case, which was discussed already in 
Ref. 

For the square lattice in D = 2 the values < k < y/2/2 are allowed. However 
if k > h only that part of a circle with radius k contributes which lies inside the 
square Brillouin zone, as shown in the inset of Figure |l|. We obtain 

f2wk iffee[Q;§] 
S2( )_ }27rfc-8fcarccos^ iffce[±;^]' ( } 

and by integrating with respect to k, 

m f^ 2 iffce[0;i] 

V2[> \^ 2 + V4F^T-4Parccos^ if ke [\;% ' 1 ' 

These functions are shown in Figure El 

In D = 3 the Brillouin zone of the simple cubic lattice allows the values < 
k < \/3/2, and more cases must be distinguished in the calculation of S3(fc). We 
omit geometric considerations and instead employ the recursion formula for S£>(fc) 
on the hypercubic lattice, 

tniii(fc,*p) 
max(0,fc 2 -i)^ 

which can be obtained in a straightforward fashion from (Q). Applying it to the 
above expression for S2(fc) we obtain, after a lengthy calculation, 



s 3 (k) = < 



ink 2 ifke [0; |] 

6nk - 8nk 2 if k G [\; % 

6irk - 87rfc 2 + 48k 2 arctan Jl--^jp 
- 24fc arctan ^/Ak 2 - 2 if k e ; & 



(15) 
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Fig. 1. Surface S2(k) and volume V2(k) of a circular Fermi sea with Fermi vector k on the square 
lattice. Inset: If the Fermi surface has radius k < h it lies completely within the first Brillouin 
zone. If k > jr the inside part contributes 2(pk in each quadrant to S2(k), where cost/) = l/(2fc), 
leading to Q. 



and one more integration yields 



37Tfc 2 - 



§**» 



if fc 6 [0; |] 
if fee [I; ^] 



1 

2A 3 



(16) 



3nk 2 - §7rfc 3 - f + 16A: 3 arctan yT 

+ \/Ak 2 - 2 - (12fc 2 - 1) arctan V4fc 2 - 2 if fee p^;^] 

These functions are displayed in Figure |2[ 

As a simple example, let us consider the rectangular (constant) density of states 

N(u) = l, 0<lj<1, (17) 

where n{n) = n, and hence e(fc) = vrj(k). Thus the plot of vo(k) in Figures 0and@ 
also represents the shape of the dispersion that will yield the density of states (fL7h . 
Next we apply the method to the model density of states (|4|). The band filling is 
now obtained as 

1 



= 1 - 



ira 



arccos fi 



2 \/l — a 2 arctan 



(l-a)(l- M ) 



(1 + o)(1+m) 



(18) 



which for the Bcthe lattice (a = 0) reduces to n(/x) = \ + — f-i 2 + arcsin/x]/7r. 

The resulting dispersion relations are shown in Figure [| Their curvature depends 
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Fig. 2. Surface S3(fc) and volume v${k) of a spherical Fermi sea with Fermi vector k on the simple 
cubic lattice. 



strongly on the asymmetry parameter a, and also on D. 

Another quantity of interest is the hopping amplitude between two lattice sites 
connected by a lattice vector R = (i?i , . . . Ru) with integer components. It is defined 
as (note the conventional overall minus sign) 



t(R) 



dki- dk D e( v /fcf 



(19) 



We calculate the hopping amplitude (|l9| ) for the one-dimensional dispersion (Fig- 
ure |§|a) . For small R we find that the amplitudes for hopping between nearest and 
between next-nearest neighbors, t(l) and t(2), are of opposite sign, and |t(2)/i(l)| 
is appreciable for a close to 1; for example t(l) = 0.337, t(2) = —0.080, for a = 0.95. 
For large \R\ it can be shownQ that t(R) cx (— 1) R |i?|~ 5 / 3 . While the stabilization 
of the ferromagnetic phase was attributed to large spectral weight below the Fermi 
energy in Ref. the present calculation shows that from a complementary point 
of view (Q) can be regarded approximately as a t-t' model, which has been shown 
to be favorable for ferromagnetism by other methods. El 



4. Limit of large dimensions 

Finally we turn to the limit D — > 00. To obtain a factorized expression for S£>(fc) 
for the hypercubic lattice we transform the argument of the delta function and use 
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Fig. 3. Dispersions in D = 1, 2, 3 leading to the model density of states (^) for several values of 
the asymmetry parameter a. 



the representation S(x) = dt e itx /(2ir), 

1 i 

2 2 

s D (k) = 2k Jdh-J dk D 8(kl 



(20) 



dte~ itk g (t) D , 

V 2 A np itq 2 



(21) 



where we defined g (t) = dg e w « . The maximum value for k is y/D/2; it is 

useful to introduce scaled quantities, K = 2k/yD and sd(k) = sd(k>vD /2)yD /2, 
Vd{k) = v d (k^/D/2), with the property / d«S£)(/c) = Vd(1) = 1. We change 
variables to s = tD/4, 



(22) 



This equation is valid for all D. From the expansion go (t) = l+it/12 — i 2 /160+O(i 3 ) 
we obtain for large 13 

IS S \ 



go(%) D = exp 



J(l + 0(^)). 



(23) 



3 10£>, 

which, when substituted into (|22j), leads to a sharply peaked Gaussian for large D, 

2 k 



s D (n) 



'2-KG 



exp 



|K ' J ^(1 + 0(1)), 



2a 2 



(24) 
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with the abbreviation a = 1/ 
tually approaches a delta function, Sd(k) 



5D. We find that for large D the function (re) even- 
l/\/3) and sd(&) 



<5(re-l/\/3) and s^fc) - S(k-y/D/12). 
This surprisingly simple result means that in this limit the Fermi sea is essentially 
located in a narrow hyperspherical shell (cut off at the Brillouin zone boundaries) at 
roughly half the maximum value of k. Interestingly, the Fermi surface has the same 
average radius as for the tight-binding dispersion with nearest-neighbor hopping.c3 
The Fermi surface in the latter case is distributed with a Gaussian in k - 1/V3 
instead of k 2 — 1/3; these distributions become equivalent for D —> oo. 

Next we obtain the asymptotic expression for vd(k) from (E4T) to the same order, 



v d (k) = 



crfc( 



eric(^4=± 



)-erfc(5^-) 



erfc( 



3V2cr 



if re e [0; j=] 
if«€[^;l] 



(25) 



This function ultimately approaches a step function, uo(k) ~ Q(k — l/y/3) and 
vo(k) ~ 6(fc — y/D/12); however further calculations must be performed with the 
asymptotic expressions ( pi| ) and ( p5| ) instead of the delta and step functions. In 
particular (^|) and ( p5| ) prescribe how to construct a dispersion that in the limit of 
infinite dimensions will yield any desired N(u) with finite bandwidth, for example 
(||) or (H). The functions Vd(k) and Sd(k) are shown in Figure ^. As D approaches 
infinity, the dispersion thus has the form of a smoothed-out step function, changed 
in shape by the application of fi{n). 

Let us consider some examples. As noted above the plot of Vjj{k) in Figure |J 
already shows the dispersion that yields the rectangular density of states (nj) . On 
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the other hand, for the special quadratic dispersion 

e(fc) = V5S(^-i) J (26 ) 
we use ( |24| ) and ([To]) to obtain the corresponding density of states 

N(cj) = — L cxp (-%r\ , -VbD/3 <lo< 2V5D/3 , (27) 



2tt 

to leading order in 1/D. Hence in the limit D — > 00 the dispersion ( p6| ) yields 
the Gaussian (||). This shows that the tight-binding spectrum for next-neighbor 
hopping on the hypercubic lattice may alternatively be represented by a quadratic 
dispersion in the limit D — > 00. 

It is instructive to calculate the hopping amplitude (|l^) corresponding to the 
dispersion (||) in the limit of large D. We change variables to |fe|, introduce the 
scaled dispersion e(n) = e(KyD/2), and perform a similar factorization as before, 
with the result (valid for all D) 

1 _ 00 D 
t{R) = -J dK ^ Jds e^ 2 J! -y;; (28) 

-00 3 = 1 

where we introduced gn(t) — J^ 2 dqe ltq2+27riqR . Note that by setting R — in 
fl2S| ) the center of mass of the band, i(0) = f UmBX du uj N(lu), is correctly recovered. 
Proceeding to the limit D — > 00 we first observe that gii(t) = —\t(—l) R / (2n 2 R 2 ) + 
0(t 2 ) for integer R ^ 0. The hopping amplitude along a coordinate axis is then 
calculated to leading order in 1/D as 

tiORi) = i((i?i,0,...,0)) 



f-l) fll /", kc(k) , , 







1 2 

IS I — — K 



3 / 10L> 



W)Af_!£M). (29) 





Using (|24|) the remaining integral is calculated by partial integration and yields 
y/3e'(l/y/3), which is evaluated via (g), 

e'(«) = M«) = tt. • ( 3 °) 



To leading order in 1/Z) we have sd (1/ \/3) = \/10D/3tt and I'd (1/ V&) = \ , whence 



= «5 (3D 



D R{ 7rliV( M (i)) 

The numerical factor if contains the value of the density of states at the Fermi 
energy for a half-filled band, /i(^). For fixed R\ this type of hopping takes place from 



10 M. Kollar 



each site to 2D other lattice sites, and the hopping amplitude is scaled proportional 
to the inverse square root of this coordination number, which is also the required 
scaling for a tight-binding band with hopping to nearest neighbors £1 An important 
difference is that in the present case the hopping is long ranged and slowly decaying. 
It is also apparent that t(R) has higher symmetry than the crystal lattice, which 
is a consequence of the spherical symmetry of eg. This spherical symmetry is also 
responsible for the large amplitude for long-range hopping, as any type of finite- 
range hopping would lead to a dispersion relation with only a discrete symmetry. 

In similar fashion we calculate the hopping amplitude between two sites con- 
nected by a lattice vector with two nonzero entries, £2(^1,-^2) = t((Ri, R2, 0, 

...,o)), 

HRi,R 2 )-p mm , t 2 - ^ N{ ^ ))Z ■ (32) 

The value of t\ now also involves the first derivative of the density of states; for 
example, t\ = -0.080 and t\ = 0.839 for (@) with a = 0.95, but t\ = for (§). 
The number of sites connected by hopping of type (|32|) is proportional to D 2 , the 
square root of which again appears in the scaling factor as expected. We conclude 
that although each hopping amplitude is scaled with dimension in the same way 
as in the tight-binding case, the present long-range hopping amplitudes combine 
to yield an energy band with a finite width, in contrast to the infinite tails in the 
tight-binding spectrum. 

Finally, note that the nesting symmetry, eg = — eg +( j, which results if hopping 
takes place only between different sublattices of a bipartite lattice, is absent for a 
spherically symmetric dispersion. Here k + Q is taken modulo the first Brillouin 
zone and the antiferromagnetic wave vector is Q = (|, 5, . . . , 5) in our notation. 
For example, the spherically symmetric representation ( p6| ) of the Gaussian density 
of states yields 



! = 1 



for \ki\ < |. Hence eg 7^ — £j: + q for most of the Brillouin zone, in particular on or 
near the Fermi surface, where |fcj| = ^= + 0(D^ 1 / 2 ). That the nesting condition 
is not fulfilled can also be seen from the fact that sites on the same sublattice are 
connected by the long-range hopping amplitude ( |3l| ) . In general, this frustration is 
expected to suppress antiferromagnetism in the Hubbard model, which is of interest 
in particular for the study of the Mott-Hubbard metal-insulator transition in the 
paramagnetic phase. It seems that spherically symmetric dispersion relations offer 
a conceptually and calculationally simpler way to remove the antiferromagnetic 
phase than, e.g., "fully frustrated" hopping on the Bethe latticed 
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5. Conclusion and outlook 

In conclusion we presented a method that, for any crystal lattice, maps a one- band 
density of states to a dispersion eg that only depends on \k\. This procedure un- 
physically enlarges the usual invariance of the dispersion under the point group 
of the crystal to a continuous spherical symmetry. As a consequence, the nesting 
symmetry typical for nearest-neighbor hopping on bipartite lattices is absent. Fur- 
thermore, the dispersion may have kinks if the specified density of states has van 
Hove singularities that are inappropriate for the dimension D of the representing 
lattice. On the other hand, the spherical symmetry may allow for analytic calcula- 
tions of DMFT one-particle quantities in the homogeneous phase that depend on 
the dispersion only through the density of states, such as the Green function of the 
Hubbard model in the vicinity of the Mott-Hubbard metal-insulator transition. 
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